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We show that dark matter haloes, in n-body simulations, have a boundary layer (BL) with 
precise features. In particular, it encloses all dynamically stable mass while, outside it, dy- 
namical stability is lost soon. Particles can pass through such BL, which however acts as a 
confinement barrier for dynamical properties. BL is set by evaluating kinetic and potential en- 
ergies (T(r) and W(r)) and calculating 1Z — —2T/W. Then, on BL, 1Z has a minimum which 
closely approaches a maximum of w — — d\og W/d log r. Such Rw "requirement" is consistent 
with virial equilibrium, but implies further regularities. We test the presence of a BL around 
haloes in spatially flat CDM simulations, with or without cosmological constant. We find 
that the mass M c , enclosed within the radius r c , where the Rw requirement is fulfilled, closely 
approaches the mass Mdyn, evaluated from the velocities of all particles within r c , according 
to the virial theorem. Using r c we can then determine an individual density contrast A c for 
each virialized halo, which can be compared with the "virial" density contrast A„ ~ 178 0^ 
(f2 m : matter density parameter) obtained assuming a spherically symmetric and unperturbed 
fluctuation growth. The spread in A c is wide, and cannot be neglected when global physical 
quantities related to the clusters are calculated, while the average A c is ~ 25 % smaller than 
the corresponding A„; moreover if Mdyn is defined from the radius linked to A„, we have a 
much worse fit with particle mass then starting from Rw requirement. 



1 Introduction 

Finding the volume where cluster materials are in virial equilibrium might seem a tough and 
somewhat ambiguous task, as clusters are not surrounded by vacuum and shall be however 
confined by unbound background materials. Here we show that, in general, the transition from 
cluster and external materials is well defined and takes place just outside a boundary layer 
(BL), which, by itself, is in virial equilibrium with no pressure pushes, but above which virial 
equilibrium is lost soon. Particles may pass through such BL, in both directions, but its overall 
virial balance keeps constant in time. We find such BL assuming a spherical symmetry. At 
a distance where the cluster ends up into unbound materials, such requirement might seem 
excessive. Rather surprisingly, however, we shall see that this is not a prohibitive restriction. 

Previous work on this point was based on the assumption of isolated fluctuation growth, 
starting from and preserving spherical symmetry (Gott & Rees 1975, Lahav et. al. 1991, Eke 



et al. 1996, 1998, Brian & Norman 1998). Results obtained in this way were then applied to 
simulation analysis. On the contrary, in our work, the assumption of isolated growth is absent, 
while spherical symmetry is tested, a posteriori, on simulated clusters. In this way we find that 
actual clusters closely approach properties predicted under such assumption. 



2 How to find BL 



Let us define integral potential and kinetic energy and their virial ratio: 
2T« r) = £ mv? , W« r) = -£ ^ , K« r) = . (2.1) 

Similar quantities can be defined for r-intervals. The r dependence of W can be outlined by 
performing the sums in eq. (2.1) starting from outside: 

W(<r) = -E E ^ = "E *(»•*)■ ( 2 - 2 ) 

i (n<r) 3 ( rj <n) %3 i (ri<r) 

Here the last term is the very definition of Z. Quite in general, a volume integral of p{r) 
increases with r. Henceforth, Z(r) = C(r/f)~ w , where C is a normalization constant evaluated 
at an arbitrary position f and w(< 1) shall however depend on r. Accordingly, it must be 

r^-(r) = -[w + rw' In Cz)]Z(r) . (2.4) 
dr r 

Let us now suppose that there exists an interval Ar = r + — r_, (i) in virial equilibrium, (ii) 
subject to no pressure pushes inwards or outwards, (iii) inside which w is constant. Owing to 
the conditions (i) and (ii), in this layer 

E mvf - r~( ri ) = . (2.5) 

i (r^eAr) 

Taking then into account the extra condition (iii), i.e. that w is constant, we can use eq. (2.4) 
with w' = and the expression for the virial ratio in A r to obtain that 

K = w , — = (2.6) 
dr 

all along the interval Ar. Viceversa, if the eq.s (2.6) are both fulfilled in a layer of depth Ar, 
such layer is at rest and in virial equilibrium. We thus define boundary layer (BL) a region of 
depth Ar satisfying eq.s (2.6). It can also be shown that no further materials can be in virial 
equilibrium if we required that w is maximum in Ar and therefore w' < for r > r + . 



3 Application to numerical simulations 



We seek BL for clusters, in simulations. We use two main sets of simulations, performed with 
different codes. The former set comprises large volume simulations (box side=360/i~ 1 Mpc), run 
with a parallel AP3M code developed by Gardini et al. (1999). They describe a tilted CDM 
(TCDM) model (n = 0.8) and a ACDM model (Q m = 0.35, ft A = 0.65) and were also used 
in Maccio et al. (2001; simulations A and B respectively). This set of simulations is meant to 
provide a significant statistics, although their mass (180 3 particles) and force resolution (~ 40kpc 
P.E.) is limited. The second set of simulations aims to provide high-resolution clusters, using 



the PM ART code (Kravtshov, Klypin & Khokhlov, 1997) and/or the public parallel tree-code 
GADGET (Springel et. al 2001). Here we resolve clusters with more than 300.000 particles 
within a radius of 2 h' 1 Mpc, with a mass resolution of 1.2 x lQ w £l m h~ l M & . With these codes 
we performed simulations C and D of ACDM (fi m = 0.3 , ft A = 0.7) and TCDM (n = 0.8) 
models, obtaining 6 high resolution clusters for each model. 

To set the position of BL, in each dark matter halo, we first calculate the virial ratio 1Z in 
successive layers; then, we select points where 1Z has a minimum and (nearly) intersects w (so 
to fulfill eq.s 2.6). Such coincidence is shown in Fig. 1 for the best and worst cases we found in 
the simulations C and D. Even in worst case, the coincidence between a minimum of 1Z and a 
maximum of w is neat. In simulations A and B, w is not so well traced, as it is obtained from 
differentiation, and is only used to select among minima of 1Z, when we have several of them. 
The procedure however works in 97% of cases, allowing to find a precise BL. 

Once the sphere confining cluster materials is set, we can evaluate the density contrast 
A c and the mass M c thereinside. In Fig.s 2 and 3, points give A c and M c for all clusters in 
simulations A and B, respectively. They show that the spread of A c values is fairly wide. By 
subdividing the M c abscissa in intervals of constant logarithmic width, we evaluate the average 
density contrast in each of them, to seek systematic trends with mass. 

Owing to the spread of A c values, their averages are still subject to a significant uncertainty, 
shown, at the 1-a level, in the plots. There seems to be no evidence of any peculiar trend 
of density contrasts with mass apart, perhaps, a modest indication of an increasing density 
contrast, at very high scales, in ACDM. It is therefore licit to consider the overall average 
among A c 's. Such average is indicated by the continuous horizontal line and compared with the 
"virial" density contrast A„, as given by eq. (1.1). In both cases, A^ (dotted line) is well inside 
the range of the density contrasts we found; the average A c , however, in both cases, is smaller 
than A„ by ~ 25 %. 

A critical result of our analysis is however shown in Fig.s 4 and 5 (for simulations A and 
B, respectively). In simulations, cluster masses can be comfortably evaluated by summing up 

particle masses and we compared this mass extimate with M^ yn = ^'g iv ' c , evaluated 

averaging over the velocities of the particles within r V:C (r v is the radius whereinside the density 
contrast is A„). This comparison yields a reasonable coincidence between M v and Md yn as is 
shown in Fig.s 4 and 5 by the dashed histograms (which confirm the slight excess of M^yn vs - 
M v already noticed by previous authors). In the same plots we also report the much better fit 
between M^ yn and M c , obtained on the basis of the setting of the BL (continuous line). The 
average values of M^ n /M c are ~ 0.97 ± 0.03 for both models. 



4 Conclusions 

Our results are based on the mathematical definition of a boundary layer (BL), bordering a region 
satisfying the virial theorem, although its density contrast has different values in different haloes. 
The BL is found in numerical simulation, around each cluster, and it is shown to separate bound 
cluster materials from the surrounding medium. Boundary layers can be found in simulations 
with sufficient dynamical range, and have been confirmed by inspections of higher resolution 
simulations. The spread of the density contrasts of individual clusters has been found to be 
fairly wide and, in average, A c are smaller by ~ 25% when compared to the so-called "virial" 
density contrast (A„ ~ 178 f^' 45 ). No evidence of a systematic trend of A c with M c was found. 
The fit between dynamical masses and particle masses, obtained with our technique, is however 
far better than the one obtained starting from a fixed virial density contrast. 
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